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Abstract — Period doubling bifurcation in buck 
converters is studied by using the harmonic bal- 
ance method. A simple dynamic model of a 
buck converter in continuous conduction mode 
under voltage mode or current mode control is 
derived. This model consists of the feedback 
connection of a linear system and a nonlinear 
one. An exact harmonic balance analysis is 
used to obtain a necessary and sufficient condi- 
tion for a period doubling bifurcation to occur. 
If such a bifurcation occurs, the analysis also 
provides information on its exact location. Us- 
ing the condition for bifurcation, a feedforward 
control is designed to eliminate the period dou- 
bling bifurcation. This results in a wider range 
of allowed source voltage, and also in improved 
line regulation. 

Keywords — bifurcation, DC-DC converter, 
harmonic balance, instability 

I. INTRODUCTION 

Several authors have investigated the occurrence of pe- 
riod doubling bifurcation in DC-DC converters (Deane 
and Hamill, 1990; Hamill et ai, 1992; Tse, 1994; Fos- 
sas and Olivar, 1996). A period doubling bifurcation 
entails loss of stability of the nominal operating con- 
dition, and as such is undesirable. Moreover, a pe- 
riod doubling route to chaos could be signaled by such 
a bifurcation, eroding the performance of the circuit. 
This has led to preliminary investigations of methods 
for prevention of period doubling bifurcations in DC- 
DC converters for a single value of source voltage in 
(Podder et al, 1995; Poddar et ai, 1998). 

In this work, analysis and control of period dou- 
bling bifurcation of a buck converter in continuous- 
conduction mode are considered. A continuous-time 
feedback system model is used as the basis for a har- 
monic balance analysis of period doubling bifurcation 
as well as control design for preventing the onset of 
period doubling. The model separates the nonlinear 



switching action of the converter from the linear fil- 
tering action. Because the model resolves dynamics 
within switching intervals, it is more accurate than the 
traditional averaged models (Middlebrook and Cuk, 
1976). 

In addition, the model proposed here is valid both 
for current mode control and voltage mode control, 
and can be applied to larger range of source voltage. 

An exact harmonic balance analysis is used to ob- 
tain a necessary and sufficient condition for a period 
doubling bifurcation to occur. If such a bifurcation 
occurs, the analysis also provides information on its 
exact location. Using the condition for bifurcation, a 
feedforward control is designed to eliminate the period 
doubling bifurcation. This results in a larger range of 
allowed source voltage, and also in improved line reg- 
ulation. (Line regulation entails regulating the output 
voltage close to a constant for different values of source 
voltage.) 

Harmonic balance analysis of period doubling bifur- 
cations has been pursued for general nonlinear systems 
in (Genesio and Tesi, 1992; Piccardi, 1994; Tesi et al, 
1996). Unlike the approximate but general treatments 
given in these references, in this paper an exact har- 
monic balance analysis is performed. This is made 
possible by the special structure of the buck converter 
model. For other converters, a more typical approxi- 
mate harmonic balance analysis can be performed. 

It is noteworthy that first-order harmonic balance 
analysis has previously been applied in power electron- 
ics in the context of small-signal modeling (Yang, 1994; 
Groves Jr., 1995). 

The remainder of the paper is organized as follows. 
In Section II, a simple dynamic model is proposed for 
the buck converter in continuous conduction mode un- 
der voltage mode or current mode control. In Section 
III, the harmonic balance method is used to study pe- 
riod doubling bifurcation for the buck converter. In 
Section IV, feedforward control is used to eliminate 
period doubling bifurcation and improve line regula- 
tion. In Section V, an illustrative example is given. 
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Figure 1: Buck converter under voltage mode control 
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Figure 2: Dynamic model for voltage mode control 



Conclusions are collected in Section VI. 

II. SIMPLE DYNAMIC MODEL 

A simple dynamic model is proposed for the buck con- 
verter in continuous conduction mode under voltage 
mode control or current mode control. 

A. Voltage Mode Control 

A buck converter under voltage mode control is shown 
in Fig. 1. The source voltage and reference voltage 
are assumed constant, and are denoted by V s and V r , 
respectively. The output voltage is v Q . The voltage 
across the diode is denoted by Vd- Output signal of the 
error amplifier is denoted by y(t). The ramp signal, 
denoted by h(t) = V t + (V h - Vi){± mod 1), is a T- 
periodic function (with h(0) = Vi and h(T) — Vh). 
The switching period is T, the switching frequency 
is f s = 1/T, and the angular switching frequency is 
uj s = 2nf s . 

The switching operation in continuous conduction 
mode (with leading-edge modulation) is as follows. 
When y(t) < h(t), the switch is on, the diode is off, and 
v<i = V s . When y(t) > h(t), the switch is off, the diode 
is on, and Vd = 0. Thus the signal Vd(t) is a square 
wave, and the switch/diode combination can be mod- 
eled as a controlled square wave generator. In this 
paper, leading-edge modulation is assumed. Similar 
analysis can be readily applied to the case of trailing- 
edge modulation, where the switch is on for y(t) > h(t) 
and the switch is off for y(t) < h(t). 

The output filter (L and C with equivalent series 
resistance (ESR) R c ) and the load (R) in Fig. 1 form 
a low-pass filter with transfer function 



V ' LC(1 + ^)s 2 + (| + R c C)s + 1 

For a buck converter with a more complex power 
stage (for example, with a second output filter), the 
model remains valid, but with a more complex transfer 
function G\{s). 

Generally the error amplifier is linear and is driven 
by the signals V r and v a . Its output can be represented 



as 

y(t) = gV r + (g 2 *v )(t) (2) 

where g is a gain constant, gi (i) is an impulse re- 
sponse function and * denotes convolution. Denote 
the transfer function associated with gi($) by Gi{s). 
Since Gi(s) depends on the control scheme, no further 
restriction is placed on it. 

Thus the buck converter in continuous conduction 
mode can be modeled as a system block diagram shown 
in Fig. 2. The output of the controlled square wave 
generator is V s if (h(t) — y(t)) is greater than 0, the 
output is if (h(t) — y(t)) is less than 0. 

B. Current Mode Control 

For a buck converter under current mode control as 
depicted in Fig. 3, the switch operation is different 
from that in voltage mode control. In current mode 
control, the switch turns on at each clock pulse and 
turns off at instants when y(t) = h(t). Here h(t) is 
a slope-compensating ramp. The switch/diode combi- 
nation can also be modeled as a controlled square wave 
generator. Therefore, a buck converter under current 
mode control can be modeled by the system block di- 
agram in Fig. 4, where a linear transfer function 

GM = 

Vd{S) 

R S (RC(R + R c )s + l) 
RLC(R + R c )s 2 + (L + RR c C)s + R 

is placed in the current feedback path. The remaining 
notation is the same as for the case of voltage mode 
control. 

C. Unified Dynamic Model 

Despite the differences in switching operation in volt- 
age mode control and current mode control, they can 
be modeled in a unified setting. Both Fig. 2 and 
Fig. 4 can be further simplified as Fig. 5, where G(s) = 
Gi(s)G 2 (s) for voltage mode control and G(s) — Gi(s) 
G 2 (s) — Gi(s) for current mode control. For either 
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Figure 3: Buck converter under current mode control 
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Figure 5: Unified dynamic model of buck converter 
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Figure 4: Dynamic model for current mode control 



case, the controlled buck converter is a combination 
of a linear system G(s) and a nonlinear square wave 
generator. 

In the remainder of the paper, more emphasis is put 
on voltage mode control. The results extend readily 
to the case of current mode control. 

III. DETERMINATION OF PERIOD 
DOUBLING BIFURCATION POINT 

In practice, the nominal operating condition is a T- 
periodic solution (in period-one mode). Representa- 
tive waveforms for y(t) and Ud(t) are shown in Fig. 6. 
At the switching instant t = d, 



y(d) = h(d) 



(3) 



When period doubling bifurcation occurs, a 2T-periodic 
solution arises from the original T-periodic solution. 
Representative 2T-periodic waveforms for y(t) and Vd(t) 
are depicted in Fig. 7. Switchings occur at t = d — 5 
and at T + d + 5, where S is a small parameter that 
vanishes at the bifurcation point. From the switching 
conditions at these two instants, it follows that 



y(d-S) = 
y(T + d + S) = 



h(d-S) 
h(T + d + S) 



(4) 
(5) 
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Figure 6: y(t), h(t) and Vd{t) in period-one mode 



The harmonic balance method is a tool that can 
be used to analyze periodic solutions in nonlinear sys- 
tems. In the buck converter, the nonlinearity results 
from the switch. In the steady-state, Vd is a periodic 
signal and can be represented by a Fourier series. By 
"balancing" the equations above (written in Fourier 
series form) at the switching instants, a condition for 
existence of a periodic solution is derived. 

The harmonic balance method is applied here to de- 
termine conditions for a period doubling bifurcation 
to occur. The basic idea is that at such a bifurca- 
tion point (shown in Fig. 8), a period-one mode and 
a period-two mode coalesce. By invoking conditions 
for the existence of each mode, the period doubling 
bifurcation point can be determined. In the following, 
the source voltage V s is used as the bifurcation param- 
eter. The critical value of the bifurcation parameter, 
denoted by V Sj *, is determined. If a different parame- 
ter is taken as the bifurcation parameter, the approach 
is similar. 

A. Harmonic Balance of Period-One Mode 

In the period-one mode, v<i(t) is periodic with angu- 
lar frequency co s , and can be expressed as the Fourier 
series 



where 



w d (t)= c ^ ejnUat 

n— — oc 
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(6) 



j2mr 



Let the duty cycle be D c . Then D c = 1 - d/T. The 
average value of Vd(t), denoted as [i>d(i)]AVE, is 

M*)]ave = co = (1 - d/T)V s = D C V S (7) 
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Figure 7: y(t), h(t) and i>d(i) in period- two mode 
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Figure 8: Period doubling bifurcation 



This agrees with the standard result from the aver- 
aging method. 

From Fig. 2, the signal at the output of the error 
amplifier is 

y (t) = gV r + (g 2 *v )(t) 

OO 

= gV r + Cne^GijnWs) (8) 

n— — oo 

Using Eqns. (8) and (3), V s can be written in terms 
of d as 



V. 



h(d) - gV r 



(1 - |)G(0) + (i)ImE~ ! i==G(jnu; 



(9) 



B. Harmonic Balance of Period- Two Mode 

Similarly, in the period- two mode, y(t) is 2T-periodic 
and can be represented as the Fourier series 

oo 

y(t)=gV r + £ c „ e ^G(^i) (10) 



where 



(^)e-^sin(2fi),ifnis odd 



Subtracting (4) from (5) and substituting (10) for y 
gives 

,V h -Vi 

S —I— = 

V °° 1 1 

( V )RC( " £ ^k~i G ^ k 2 H) Sin((2fc ~ 1)Ws<5) 
fc— l 

oo 

+ V - r G(jfcw s )(sin(2fcw ;s( 5) - 2e^ d sin(fcw s <5)) 
fc=i 2k 

Solving for V s gives another expression for V s in 
terms of d and (5. See (Fang, 1997) for the detailed 
form. 

C. Determination of Bifurcation Point 

Recall that 6 is small and that, if a period doubling 
bifurcation occurs, then S = at the bifurcation point. 
The critical value of V s at the bifurcation is then de- 
termined: 



v.. 



Vh-Vi 
2 



m- rr=i G u( k - jH) + a - e^^cc/^)] 



(ii) 

The critical values V^* and rf* can be obtained graph- 
ically by plotting Eqns. (9) and (11) on the same axes. 
The intersection (V^,*, d*) of these graphs (if it oc- 
curs) is the period doubling bifurcation point. In- 
deed, a necessary and sufficient condition for a pe- 
riod doubling bifurcation to occur is that these graphs 
intersect. (Period doubling bifurcation occurs when 
the source votage reaches V s ^, or duty cycle reaches 
d„/T.) 

The denominator of (11) can be approximated by 
the term that involves G with the smallest argument, 
namely, Re[G(ju; s ) — G(^)]. So a good estimate for 
the critical value V s * is 



Vh-Vi 
2 



Re[G{ju s ) - G(^)] 



(12) 



Further bifurcation to the period-four mode can be 
analyzed similarly and is obmitted here. Two special 
cases are considered next. 

Example 1 Consider a buck converter under voltage 
mode control with G2(s) = gi and g = —gi, shown in 
Fig. 9. Under the conditions \j\f~LC < ui s , R c < R, 
and R C C <C 1 (which are generally true), the critical 
source voltage (Eqn. (12)) can be further simplified as 



V s . 



V h 



6g 



Rc 
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)LCw 2 a 



(13) 



(^)(e" i ^ cos(^) - e^), if n is even 



It shows that larger values of uj s , L, C, and ramp 
amplitude (Vh — Vi) lead to larger stable range of source 
voltage, while larger feedback gain gi leads to smaller 
range of source voltage. Also ESR R c (compared to 
R) is generally small and has little effect on stability. 
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Figure 9: System diagram for Example 1 



Critical values of other parameters can be obtained 
similarly from Eqn. (13). For example, the critical 
feedback gain gi t * is close to 



6K 



R 



(14) 



Example 2 Consider a buck converter under cur- 
rent mode control with open voltage loop (G*2(s) = 
0). Under the conditions u s >• (l/^/LC, R c /L, and 
1/RC), the critical source voltage (Eqn. (12)) can be 
further simplified as 



V. 



Rr. 



6R S 



(15) 



With slope compensation and larger values of uj s 
and L, the stable range of source voltage is larger. 
However, different from voltage mode control, larger 
ESR R c leads to smaller range of source voltage. 

These two examples shows that if a control scheme 
is designed appropriately to adjust the amplitude of 
the ramp, the period doubling bifurcation can be pre- 
vented. This is pursued next. 

IV. FEEDFORWARD CONTROL 

Common objectives of controller design for a DC-DC 
converter are stability, fast transient dynamics, line 
regulation, and load regulation. In some demanding 
situations, not all of these objectives can be achieved 
simultaneously due to nonlinear nature of the con- 
verter. In this section, feedforward control is used to 
prevent period doubling bifurcation (instability) and 
to achieve line regulation. Transient dynamics and 
load regulation can be improved by redesigning the 
voltage feedback loop, however, this is not the focus of 
this paper. 

Here a feedforward control scheme is proposed to 
adjust the ramp signal h(t) by setting Vi = kiV s and 
Vh = khV s , where fc; and kh are the feedforward gains. 
In this control scheme, the amplitude of the ramp, 
Vh — Vi = (kh — ki)V s , is proportional to the source 
voltage. 



A. Prevention of Period Doubling Bifurcation 

Define the following function of d: 



1, 



fe=i 



H(d) =: 2Re[- VGOXfc-^V s )+(l-e^ d )G(jfc 



(16) 

Denote its maximum by H max and its minimum by 
The condition for bifurcation in Eqn. (11) becomes 



H{d) = k h - h 



(17) 



If the values of kh and k\ are chosen to satisfy (kh — 
h) > i? max or (kh - k) < H min , the period doubling 
bifurcation is prevented because the bifurcation con- 
dition is never met. 

B. Line Regulation 

In the last subsection, one has two degrees of freedom 
in choosing the values of kh and fc;. This freedom can 
be used to achieve another objective besides prevent- 
ing period doubling bifurcation. A common objective 
in DC-DC conversion is line regulation, which is ad- 
dressed in this subsection. 

First, an equation related to the duty cycle D c = 
1 — d/T is derived. Assume the switching frequency 
is high enough that the high order terms in Eqn. (9) 
can be ignored because of the low-pass nature of G(s). 
Then Eqn. (9) becomes 

Vi + (V h -Vi)±-gVr 



V 



(1-|)G(0) 
Solving this equation for d/T gives 



d 
T 



G(0)V s -V+ gV r 



(18) 



(19) 



G(0)V s -V + V h 

Next, the average output voltage, denoted as [w ]ave, 
can be related to the source voltage V s as 

[w ]ave = Ki]ave, irom Fi g- 2 and G i(°) = 1 

= (1-|)K, from Eqn. (7) 



= ( 



V h - gV r 



G(0)V s -Vt + Vh 
k h V s - gV r 



)V S , from Eqn. (19) 

G(0) - h + k h (20) 

Two approaches make output voltage independent 
of source voltage: large G(0) (such as integral control) 
or simply setting kh = 0. The second approach is 
pursued here. From Eqn. (20) with kh = 0, one has 

gv r 



h = G(0) + 



Wo] AVE 



(21) 



Thus, given a buck converter power stage (G\(s)) with 
desired average output [w q ]ave, and an error ampli- 
fier (g and G2(s)), line regulation can be achieved by 
adding a feedforward loop to adjust the ramp h(t) by 
setting Vh = and Vi = kiV s with ki as in Eqn. (21). 
For the case in Example 1, fc; — gi(l — V r /[v ] AVE ). 
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Figure 10: Bifurcation diagram of the circuit in Fig. 9 



C. Combined Period Doubling Prevention 
and Line Regulation 

If kh = and the value of ki given by Eqn. (21) satisfies 
ki < — i? max or ki > — -ffmin, both prevention of period 
doubling bifurcation and line regulation are achieved. 

V. ILLUSTRATIVE EXAMPLE 

Consider the buck converter in Example 1 with the 
same system parameters as (Hamill et al., 1992): T = 
400 fis, L = 20 mH, C = 47 fjF, R = 22 fl, V r = 11.3 
V, 9l = 8.4, Vi = 3.8 V, V h = 8.2 V, (then h(t) = 
3.8 + 4.4[y mod 1]). Let V s be the bifurcation pa- 
rameter. Here the focus is on determination of the 
bifurcation point and controller design, although the 
switching frequency is low. 

It has been shown that a period doubling bifurca- 
tion occurs at V s — 24.5 through simulation (Hamill 
et al, 1992) and through calculating the locus of the 
eigenvalues of a discrete-time model (Fossas and Oli- 
var, 1996; Fang, 1997). The bifurcation diagram is 
shown in Fig. 10. From the diagram, the stable range 
of V s is short (16 V to 24.5 V). In this range of V s , the 
output voltage varies from 11.9 V to 12.03 V as the 
source voltage varies. 

The period doubling bifurcation point can be deter- 
mined by plotting Eqns. (9) and (11) together (see Fig. 
11). The intersection (V s ,*, d*) =(24.5, 2.04 x 10~ 4 ) of 
these graphs is the period doubling bifurcation point. 

To show the effect of ESR R c and T on the critical 
source voltage V s ,*, another two conditions (R c =l fi 
and T=250 fis) are considered. The exact and esti- 
mated V Sy * are shown in Table 1. It shows that larger 
uj s leads to larger V S} * and ESR R c has little effect on 
V s ,*. It also shows that Eqn. (13), expressed explicitly 
in terms of system parameters, gives close estimate of 

Vs.*. 

In the following, a feedforward scheme is designed 
and added to the original controlled buck converter to 
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Figure 11: Dashed line: Eqn. (9); solid line: Eqn. (11) 



Table 1: Determination of V Si * for different T and R c 



Switching period T (fis) 


400 


400 


250 


ESR R c (O) 





1 





Exact V St * by Eqns. (9), (11) 


24.5 


24.9 


49.5 


Estimate of V St * by Eqn. (12) 


20.2 


22.4 


51.8 


Estimate of V s ,* by Eqn. (13) 


20.2 


21.2 


51.8 



achieve two objectives: (i) a wider V s operating range 
(16 V to 35 V) without period doubling bifurcation 
and (ii) line regulation with [i> ]ave = 10 V. 

A plot of H{d) is shown in Fig. 12 (here, i? max = 
0.358 and if mi „ = 0.1792). Set k h = 0. From Eqn. (21), 
ki = —1.092. The condition fc; < — iJ max is satisfied. 
Thus both prevention of period doubling bifurcation 
and line regulation can be achieved. 

A feedforward control to adjust the ramp amplitude 
is designed. The ramp signal has Vi = and Vh = 
— 1.092 V s . The resulting bifurcation diagram of the 
buck converter with the feedforward control is shown 
in Fig. 13. Compared with Fig. 10, Fig. 13 shows a 
wider operating range for V s (16 V to 35 V) and good 
line regulation, with [u d ]ave = 10 V. 

Take V s = 28 V for example. The output voltage re- 
sponse during start-up (starting from (il, v a ) — (0,0)) 
is shown in Fig. 14. The output voltage is regulated 
to 10 V. The switching operation depends on intersec- 
tion of the signals h(t) and y(t). These are shown in 
Fig. 15. Different from the traditional positive ramp, 
the signal h(t) is negative. It is because the regulated 
output voltage is smaller than V r , and a negative h(t) 
is needed to compare with negative y = gi(v a — V r ). 

VI. CONCLUSIONS 

Exact harmonic balance analysis is applied to study 
period doubling bifurcation in the buck converter in 



Figure 12: Plot of function H(d) 
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Figure 13: Bifurcation diagram (feedforward control) 



continuous conduction mode. A simple and unified dy- 
namic model of the buck converter under voltage mode 
or current mode control is derived. The model is gen- 
eral that it can be applied to the buck converter with 
various configurations, such as with a second output 
filter, with ESR R c , or with a high-order error ampli- 
fier. This model consists of the feedback connection of 
a linear system and a nonlinear one. An exact condi- 
tion for period doubling bifurcation is given in terms 
of solving a pair of algebraic equations. The critical 
condition can be approximated explicitly in terms of 
system parameters, as in Eqns. (13)- (15). These con- 
ditions for buck converters show that larger values of 
switching frequency and ramp amplitude lead to larger 
stable range of source voltage. ESR R c does not affect 
stability in voltage mode control, but it does in current 
mode control. To prevent the period doubling bifurca- 
tion, a ramp-adjusting feedforward control is designed. 



Figure 14: Output response with feedforward control 
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Figure 15: h(t) and y(t) with feedforward control 



A wider operating range of source voltage is achieved, 
along with line regulation. Simulations are given to 
illustrate the effectiveness of the design technique. 
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